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Abstract 


The link-calculus has been recently proposed as a process calcu- 
lus for representing interactions that are open (i.e. that the number 
of processes may vary), and multiparty (i.e. that may involve more 
than two processes). Here, we apply the link-calculus for expressing, 
possibly hierarchical and non dyadic, biological interactions. In par- 
ticular, we provide a natural encoding of Cardelli’s Brane calculus, a 
compartment-based calculus, introduced to model the behaviour of 
nested membranes. Notably, the link-calculus is flat, but we can 
model membranes just as special processes taking part in the biological 
reaction. Moreover, we give evidence that the link-calculus allows one 
to directly model biological phenomena at the more appropriate level 
of abstraction. 


Keywords: Membrane Interactions, link-calculus, Brane Calculus 
1 Introduction 
In the last years, many process algebras have been adapted for modelling and 


analysing biological systems, by exploiting the similarity between distributed 
and concurrent systems and complex biological systems. Species can therefore 
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be abstractly seen as interacting processes and reactions as actions, thus 
helping the modeller to focus on the relationships among the species. As 
a dowry, process algebras provide a formal setting for describing biological 
phenomena and ready-to-use analysis methodologies and tools. 


The symbolic abstraction offered by process algebras is crucial for 
managing the complexity of biological mechanisms, especially if it can be 
joined to a flexible modelling strategy able to switch between different levels 
of abstraction. 


We focus here on process calculi that model biological compartments 
or membranes and their interactions, and in particular on Brane Calculi 
(presented in Section 5), a family of process calculi, introduced by Cardelli 
in [8], endowed with dynamically nested membranes, based on location 
mobility and awareness. Membrane interactions are described into terms 
of synchronisations between pairs of co-actions. Nevertheless, on a closer 
view, these apparently dyadic interactions are inherently multiparty. Each 
interaction not only involves the two processes triggering the interaction, by 
offering the required co-actions, but it also involves one or more membranes, 
as supporting actors. Each interaction leads indeed to a new membrane 
hierarchy. 


To better capture this multiparty nature of interactions, we resort to a 
flat calculus, the link-calculus [4] (recalled in Section 3), whose primitives 
are designed for describing and coordinating interactions that are multiparty, 
i.e. that may involve more than two processes and are open, i.e. the number 
of involved processes is not fixed or known a priori by the other partners. The 
link-calculus is an extension of the 7-calculus, where dyadic interactions 
are replaced by a more general synchronisation algebra. 


The main difference, with respect to the a-calculus, is that in the 
lLink-calculus communication actions are given in terms of links, that record 
the source and the target ends of each hop of interaction. Links can be 
indeed combined in link chains that route information across processes from 
a source to a destination. Intuitively, an interaction can be imagined like the 
composition of a one-dimensional jigsaw puzzle: it requires different pieces to 
fit together, where the pieces are the links that, together with the objects of 
communication, are offered by each involved action. Links can be joined in 
a link chain if they are to some extent “complementary”, i.e. if each process 
contributes with links that are not specified by others. According to the 
puzzle analogy, different parts of the chain can be composed separately, and, 
afterwards, assembled without overlays. 
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The link-calculus has been introduced and exploited in [4] to provide an 
encoding of Mobile Ambients [9], from which many bio-inspired calculi derive, 
included the Brane Calculi. Mobile Ambients and its descendant calculi are 
all based on compartmentalisation and on location-awareness. Ambients 
are bounded locations representing the mobile computational units in which 
agents interact. The link-calculus can naturally capture the spatial aspects 
due to ambient nesting and the multiparty essence of ambient interactions. 


Supported by the “genetic” resemblance with Mobile Ambients, we 
provide an encoding of Brane Calculi in the link-calculus, by following the 
line of [4]. Also in the case of Brane Calculi, by using the link-calculus we 
can easily handle locality, without introducing any specific operator, just 
encoding any membrane compartment as a separate process. Furthermore, 
multiparty interactions are naturally rendered. 


One of the peculiar features of Brane calculi is the design of different 
sets of primitives to capture different biological mechanisms, and to choose 
different levels of granularity for the studied biological compartments. We 
first encode the MBD version of Brane calculi, and then we introduce similar 
encodings for the PEP version of the calculus, and for a subset of the calculus 
able to handle the interaction between molecules. Our flat calculus easily 
accommodates each new set of primitives. Under this regard, our calculus can 
be seen as a sort of universal low-level language for the family of membrane- 
based calculi, offering a basic set of interactions able to simply implement 
the several more abstract primitives of Brane. This shows the expressivity 
of link-calculus, both from the linguistic and biological perspective. 


At the same time, we show that the link-calculus can also directly 
represent high-level molecular interactions and biological phenomena, by 
using exactly the same principles. Actually, an expressivity problem arises 
when we consider that, in general, biological interactions, such as reactions, 
can involve more than two elements at the time, but are modelled, by 
usual biologically-oriented process calculi, as suitable sequences of binary 
reactions. With the link-calculus, we do not have this limitation in the 
arity of participants in a reaction any longer. We can indeed — as shown in 
the next Section — directly model a biological process involving multiparty 
interactions in the link-calculus, without passing from another process 
algebra. One advantage is a tight operational correspondence between 
biological reactions and process reactions: while biological reductions are 
often encoded as sequences of interactions in ordinary process calculi, they 
can be represented as atomic interactions in the link-calculus. Furthermore, 
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in modelling we can choose the abstraction level, depending on the current 
interest. 


Related Works In [19, 20], there are other proposals for changing the 
abstraction level, using the same formalism: the probabilistic model checker 
PRISM is indeed applied to simulate biological case studies, both at the 
individual level, and at the level of interacting populations of cells. In [19] 
also membrane compartments are modelled. To allow the specification 
of multi-reactant multi-product reactions in BlenX [16], the work in [13] 
proposes an extension to model a sequence of elementary actions as if it were 
atomic. More similar in spirit to our work is instead [27], where an extension 
of the z-calculus, 7@, is introduced that provides an encoding of Brane 
Calculi. This calculus is flat, but with dyadic interactions: localisation is 
expressed by means of polyadic synchronisation, while atomic operations are 
expressed by means of constraints on the order of communications, obtained 
with priorities. 

Several approaches, among which we recall [17, 22, 24, 5, 6], present 
encodings of Mobile Ambients. The main motivation is to provide an LTS 
semantics to Mobile Ambients, as a basis for bisimulation congruence. Also, 
n [12], the authors present an encoding of the mobile ambients without 
communication into a subset of the 7-calculus, based on the separation of the 
spatial structure of mobile ambients from their operational semantics. While 
in all these works ad hoc semantics are introduced, the encoding proposed 
in [4], and the one presented here for Brane Calculi are just illustrative 
examples of application of our link-calculus. 

The basic computer interaction, where a message is sent and received 
by two different processes, is dyadic and there is a long tradition in studying 
the properties of this kind of interaction. The seminal work [26] that opens 
the study of the process algebra for modelling biological processes resorts to 
the z-calculus as it is one with a strong established theory. Until now, one 
of the main challenges in this area is how translating the classical properties 
of distributed systems into biological terms. Nevertheless, from the point of 
view of the intrinsic dynamics of biological processes, the interactions are 
not dyadic and their formalisation in terms of input-output activities is a 
strong simplification. 

Our proposal goes in the direction of defining a multiparty calculus 
suitable to compose a large variety of biological interactions, different in the 
number of their participants and on their level of abstraction (biochemical, 


A Flat Process Calculus for Nested Membrane Interactions 95 


molecular, cellular, etc.), sufficiently general to easily incorporate most of 
the relevant aspects in biology, such as locality, shapes, concentrations, etc. 


Structure of the paper. In Section 2, we give a first introduction to the 
link-calculus and we show its flexibility in modelling biological interactions, 
by giving a model of the Receptor-mediated Endocytosis mechanism. The 
example is then resumed and concluded in Section 4, after the formal 
introduction of the syntax and operational semantics of the link-calculus 
in Section 3. In Section 5, we recall the syntax and operational semantics 
of the families of Brane calculi, while in Section 6 we show how to encode 
Brane calculi in the link-calculus. In Section 7, we present the Link-calculus 
encoding of the Brane Calculus description of the Semliki Forest virus [8] as 
a concrete example of the translation. Some concluding remarks are drawn 
in Section 8. 


2 The Receptor-mediated Endocytosis Example 


Modelling complex biological phenomena is not a trivial task. Consider, 
for instance, the process known as receptor-mediated endocytosis [1], a very 
general mechanism that allows living cells to transport specific substances 
(ligands) from the external environment to the cytoplasm. 

This process, depicted in Figure 1, consists of various steps and involves 
the formation of vesicles. The first step occurs on the outer side of the cell 
membrane, and is represented by the binding between the ligands and a set 
of specialised molecules (receptors), embedded in the cell membrane. The 
formation of the ligand-receptor complex (LR) triggers further events close 
to the inner side of the cell membrane. First, one molecule of the AP-2 
complex? binds to each LR complex, thus allowing the clustering of the 
LR complexes and the subsequent link with a number of proteins called 
CLATHRIN. These molecules elicit the invagination of the part of the cell 
membrane that resides close to the receptors, and eventually leads to the 
formation of a vesicle that contains the LR complex and is coated by an 
external layer of AP-2 and clathrin molecules. Then, this vesicle moves 
along the cytoplasm towards other vesicles (endosomes). Finally, the vesicle 
coalesces with the endosome (fusion), immediately after the dissociation of 
the external AP-2-clathrin coat (uncoating). 


>The AP-2 complex consists of two large adaptins (a and $2), a medium adaptin (2), 
and a small adaptin (02). 
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Figure 1: The Receptor-mediated Endocytosis. 


We focus our attention on the first two steps: (1) the binding between 
the ligands and a set of specialised molecules (receptors) embedded in the 
cell membrane; and (2) the further events close to the inner side of the 
cell membrane triggered by the formation of the ligand-receptor complex 
(LR): (a) the binding of one molecule of the AP-2 complex (a.k.a. adaptor 
complex) to each LR complex, and (b) the following clustering of the LR 
complexes and (c), finally, the subsequent link with a number of proteins 
called CLATHRIN. 

To obtain the overall picture of this process, we need to address two 
different abstract levels, usually modelled separately: the level of membranes 
and the one of molecules. 

The flexibility of the link-calculus is exploited here to model both 
molecular and membrane interactions, without resorting to any extension 
of the language. This example will also give us the chance to smoothly 
introduce our calculus, that will be formally presented in Section 3. 

In the first part of the example, there are different possibilities for 
modelling the biological interactions, depending on the particular object of 
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study; we can choose to model the engulfing of the LIGAND: 


(a) as a single multiparty interaction, that involves the cell membrane, 
the LIGAND, the RECEPTOR, and the two molecules AP-2 and 
CLATHRIN; 


(b) or by resorting to two steps, by making explicit the binding between the 
ligands, the receptors, and the proteins AP-2, then the ligand between 
the ligand-receptor-AP-2 complex and the proteins CLATHRIN; 


(c) or by choosing a lower abstract level, not considering membranes, and 
only focussing on the underlying chemical reactions. 


Here, we follow the second solution (b), and we model the first part of 
the example as the two following link-calculus multiparty interactions: 


e the interaction among the LIGAND, the RECEPTOR, and the cell 
membrane, which is considered as an active participant; and 


e the interaction between the complex LIGAND-RECEPTOR and the 
proteins AP-2 and CLATHRIN. 


In Figure 2, there are the link-calculus specifications of the involved 
elements, whose names are written in capital letters, e.g. LIGAND or RE- 
CEPTOR, and where the subscript records the name of the compartment 
where each element lies (e.g. LIGAND}, represents the element LIGAND 
lying at top level in the membrane hierarchy). In membranes also the parent 
compartments are recorded as superscripts (as in the previous examples): 
e.g. Mer is the cellular membrane, lying at the top (of our model). 

In traditional process algebras, processes communicate by means of 
actions and co-actions that synchronise on the same channel and possibly 
exchange data. In link-calculus, interactions are multi-party and can there- 
fore involve more than two processes. Actions come with links, i.e. pairs 
“\3 that record the source and the destination ends of each hop of a com- 
munication. At the same time, actions offer a communication tuple (w), 
whose names can be used either as values (e.g. cell) or as variables (that we 
distinguish by underlining them, e.g. xj). In the cell membrane process, in 
the link part ("\...) we specify the binding elements placed on the membrane 
(in our case, the external and the internal part of the receptor), while in 
the tuple, we specify the two involved localities (top and cell). Instead, 
in the link part of the others biological elements, we specify the element 
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name, on the left side, and the element name that can be bound to, on the 
tight side (Fier ar Name, BindingElement). AS before, we specify the element 
location in the tuple. Each party offers its link and its tuple, and a multi- 
party interaction is possible when they can be suitably combined together, 
i.e. when: 


e links match two by two, e.g. ™\memRec With TT RE ot and when 


e all the parties agree on the same tuple pattern, by possibly matching 
values in corresponding position, and by binding the remaining variables 
with the corresponding values, which here represent the names of 
their own locations. For instance, in the combination of the tuples 
(COD; Dracminias Bidens Lael: ONG. (iops NEM ants Pinan, steed), whe value 
top can be assigned to the variable xj), while the value memoyt can 
be assigned to the variable tmem,,,;- 


Therefore, links are combined in a link chain, like pieces in a jigsaw 
puzzle, where each party contributes with its link. More precisely, we can 
think about S-shaped tetrominos (see the first link chain produced in our 
example in Figure 3), joined together when the labels of the edges match. 

The following first transition allows the formation of the ligand-receptor- 
AP2 complex, while the second transition produces COM PLEX ey con- 
taining the vesicle membrane and the complex ligand-receptor-AP2-clathrin. 
The definition of COMPLEX,,, is given in Figure 4. 


LIGAND op | M22? [RECEPTOR ceti| AP-2ce1t| CLATHRIN ce1i 


cell 


{ Dene Ree ie SE ae 
memRec \receptor \ap2 \T 


M*? |LIG-REC-AP2c¢1;|CLATHRIN ce1i 


cell 
bo eee Me 


M£? COMPLEX cei 


cell 


3 The Calculus of Linked Interactions 


In the previous section, we intuitively introduce the calculus of linked in- 
teractions, (link-calculus for short). Now, we formally introduce its main 
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LIGANDiop| Mich [RECEPTOR ce|AP-2ce1|CLATHRIN cet 
(whole system) 


LIGAND top = 7 \ memRec (top, Lmemours UMemin > Leell)-O 


(outside the cell) 


top A mem Rec 
Mei — recx.( \ receptor (ees MEMout,; Umemin> Leell) 


(the cell membrane) 


RECEPTOR cei = ee a te (Xtop; Tmemouts MVEMin, Leell) .LIG-REC-AP2cei 
(on the cell membrane) 


A ap2 
AP-2cetl = %P \r (Pian, Lmemouts Umemin> cell) .O 


(within the cell) 


CLATHRIN, ey & “2#r"\ , (x, cell).0 
(within the cell) 


LIG-REC-AP2ceu = *\ciathrin (cell, y).COM PLEX ceit 
(within the cell) 


Figure 2: Model of the formation of the complex 
ligand—receptor—ap2—clathrin. 


ingredients, and we briefly recall its syntax and semantics (for further details 
and examples we refer to [4]). The key feature of this language is that 
communication actions are given in terms of links. A link is a pair that 
records the source and the target ends of a communication, meaning that 
the input available at the source end can be forwarded to the target one. 
Links are combined in link chains to describe the way information can be 
routed across ends. 


Let C be the set of channels, ranged over by a,b,c, ..., and let CU{7T, * } 
be the set of actions, ranged over by a, 3,7, ..., where the symbol 7 denotes a 
silent action, and the symbol « denotes a non-specified action (i.e. a missing 
piece of the puzzle according to the analogy proposed before). 

A link is a pair = %\g: we call a the source end of £ and £ the target 
end of &. A link = \g is valid if either a, 8 ¢ * or £=*\,. In the first 
case, the link is called solid. The link *\, is called virtual. We let £L be the 
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t memRec receptor ap2 


memRec receptor ap2 t 


Figure 3: Link Composition seen as Combining Tetrominos with Matching 
Labels. 


set of valid links. Examples of non valid links are 7\, and *\, while both 
7\, and ®\, are valid links. 

A link chain is a finite sequence s = ¢1...£,, of (valid) links @; = °\g, 
such that: 


Bi, 041 €C implies B; = ai41 


1. for any i€ (L,n— 1, { BPS i yaa = 7 
i= a el 


2. if Vi € [1,n].a;,G € {7,*}, then Vi € [1,n].a; = B; = T. 


The first condition says that any two adjacent solid links must agree on their 
ends: it also imposes that 7 cannot be matched by *. The second condition 
disallows chains made of virtual links only. The empty link chain is denoted 
by «. A non-empty link chain is solid if all its links are so. For example, 
7\@\) is a solid link chain, while ™\*\°\- is not solid. 

The link-calculus processes are defined by the following grammar: 


P.Q:=0| X | &.P | P+Q| PIQ | (va)P | Pla/d] | recX.P 


where ¢ is a solid link (i.e. = \g with a, 6 # *) and t is a tuple of names. 

Roughly, processes are built over a standard syntax (with nil process 0, 
prefix ¢t.P, sum P + Q, parallel P|Q, restriction (va)P, renaming P{a/b], 
recursion recX. P, for X a process variable), but where the underlying 
synchronisation algebra is based on link chains. 

The calculus is equipped with the name passing mechanism, borrowed 
from the (polyadic) z-calculus, but slightly extended. Each link in the 
whole chain simply carries the same list of arguments, but with different 
(send/receive) capabilities. We assume channel names are admissible values, 
i.e. as in 7-calculus, we have the possibility to communicate (names of) 
means of communication. In the tuple t = (w), names can be used either 
as values or as variables. To be distinguished, variables are underlined. 
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P+0=P P+Q=Q+P (P+Q)+R=P+4+(Q+R) 

P|O=P P/Q=Q|P_ (P|Q)|R=P|(Q|R) 

(vn)0=0 (vn)(um)P = (vm)(vn)P 
recX.P = P{recX.P/P} (vn)(P|Q)=P| (vn)Q, if n¢ fn(P) 


Table 1: Congruence rules of the link-calculus. 


This mechanism allows, e.g. a form of multi-way communication, where 
all peers involved in the link chain can express arguments to be matched, 
and provide actual arguments to replace the formal ones of other peers. 
In general, the same peer can atomically (i.e. executing one single prefix) 
send some values and receive other values. For a tuple t, we let vals(t) and 
vars(t) denote the set of values and the set of variables of t, respectively. 
We say that a tuple t is ground if vars(t) = 0. For example, the processes 
P=7\,(v).P’ and Q = “\,(x).Q’ can interact by forming the link chain 
7\¢\_ with P sending the value v to Q (x will be substituted by v in Q’ after 
the interaction). As a more general example, processes P = 7\q(v,y,m).P’ 
and Q = “\,(a, w,m).Q! can interact by forming again the link chain 7\¢\, 
with P sending the value v to Q (as before), Q sending the value w to P 
(y will be substituted by w in P’ after the interaction) and both P and Q 
agreeing over the third argument m of the communication (which is matched 
in both tuples). 


As usual, (va)P binds the occurrences of a in P. Similarly to (va), 
the prefix @t.P binds the occurrences of the variables vars(t) in P. The sets 
of free and of bound names of a process P are defined in the obvious way 
and denoted, respectively, by fn(P) and bn(P). Processes are taken up to 
alpha-conversion of bound names, and we shall often omit trailing 0, and 
empty tuples, e.g. by writing “\, instead of “\,().0. In the following, given 
two set of names S' and T’,, we write S#T as a shorthand for SN T = 90. 


The operational semantics is defined in terms of an LTS, whose states 
are link-calculus processes, whose labels are pairs sg of link chains s and 
tuples g, and whose transitions are generated by the SOS rules in Table 2, 
by relying on the congruence rules in Table 1. In a label of the form sg, if s 
is solid, then g must be the empty tuple (), i.e. it is not possible to observe 
the arguments of a completed communication. 

We denote by o = [11 + U1,...,2n > Up] the substitution that replaces 
each x; with v;, and is the identity otherwise, and set vars(o) = { £1,..., Un }. 
Substitution is defined, as usual, in order to avoid name captures, and can 
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< ag. 
les — Zo t (Act) Se (Lsum) 
at.P 3 Po P+Q4P! 
P sg Pp’ P 8g Pp! 
— ) : : : — = me = vane 
(va)P ——*; (va)P’ (va)P ——, Pp 
PIX recX PIP PSP exlg)#fn(@) 
ec par 
recX. P “4 P’ PIQ = PQ 


P 4, pl o= LT 6 eee ses’ is not solid 


= (Com) 
PIQ = PQ’ 
89. ay / ex(g)#fn(Q) ses’ is solid 
a q— +Q ex(g’)#fn(P) geg’ is ground Ress 
PIQ 5 (v ex(g #9'))(P'1Q') 
89, pr = _—_ 8 _yp 

Pate abe (Ren) pep ie Q=Q (Struct) 

Pla/t] “9° P'a/o PQ’ 


Table 2: SOS Semantics of the link-calculus (rules (Rsum) and (Rpar) 
omitted). 


be applied to tuples, links, and to processes (see [4] for further details). 
We say that g is a full instance of t and write g x, t if vars(o) = vars(t) 
and g = to. 


Certain actions of the link chain can be hidden, by restricting the 
channel where they take place. Formally, for s = ¢1...€n, with 0; = \,, for 
€ [1,n], we define the restriction operation (va)s by letting® 


®This operation and the following ones are not defined in all the cases not explicitly 
mentioned. 
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(va)s = ™\(va)(3?)\..\(va)(B"_, )\bn if a1, Bn # @ 
7 foe p =a 
- o ate e 


oN 

y 

2 
bee sk 
em, 
Rene) 
ait 
|> 


Based on this definition and like in the 7-calculus, names in the tuple 
can be extruded during the communication. In the labels of transitions, we 
need to annotate positions in the tuple to distinguish between arguments 
that are taken in input (i.e. they are guessed instances), or that are extruded. 
We underline the former ones, while we overline the latter ones. A name can 
be extruded when it is not already annotated; after the extrusion, it will be 
overlined. Formally, given a link chain s, an (annotated) tuple g, (v.a)(sg) 
and (va)g are defined as follows: 


(v.a)(sg) 
(va) (wy, ..., Wn) 


(vajw = . AE Gs 8 
a ifw=a 


> |l> 
eas 
——~ 
Q 
wW 
~ 
F aii | 
——~ 

5 ON 
Q 
~S" 
SC 
~ 


We let ex(g) denote the set of extruded (i.e. overlined) names appearing in 
g, and extend the definition of substitution application to overlined names 
in the obvious way. We write a © g if the name a appears in the tuple g 
(with or without annotation). 

Two link chains can be merged if they are to some extent “complemen- 
tary”, i.e. if they have the same length, each provides links that are not 
specified in the other, and together they form a (valid) link chain. Formally, 
fore Piel On oe Clg wri ey ee end 1 \ gy for any 
i € [1,n], we define se s’ by letting: 


ses’ 
“\a @ © \ gr 


Two annotated tuples can be merged when they list exactly the same values 
in the same order, and if the values in matching positions are annotated in 
some compatible way. Formally, if g = (w1,...,Wn) and g’ = (wi,..., wi): 


(C, @ £;)...(Lne £,) 
(a0 0)\ (08) 


I> 


ae B 


I> ||> 


a iffG=x 
B ifa=x 


(ses')(geg’) 
(Wy © Wy.) Wn © Wy.) 


sges'g! 
geg’ 


A 
A 
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w fwHew =v) Vw =v) 
wew + 0- Hie =0hG. =u) Ve =on w=) 
Oo two =n (w= vAw =F) 


The SOS rules in Table 2 resemble the early semantic rules of 7-calculus. 
The main difference is that we are dealing with a multi-party form of 
communication, and, consequently, the “close” rule must be applied when 
the communication has been completed. 

Rule (Act) allows the process ¢¢.P to offer the tuple t in a communication 
on the link @. More precisely, following an early style, the actual tuple to be 
communicated must be a full instance of t (see the condition g <, t): the 
communication is that of sg, where variables appearing in t are replaced 
in g = to by actual parameters. The substitution ¢ is also applied to the 
continuation, after the transition (Pc). 

In rules (Res) and (Open), we leave implicit the side condition (va)sg # 
1. Rule (Res) is applicable whenever a ¢ g, in which case (va)g = g 4 L 
and thus (v a)(sg) = ((va)s)g. Rule (Open) models the extrusion of a. Note 
that, since (va)sg # | and a € g, we have that the only possibility for a 
to appear in g is without annotations (otherwise (va)g = L). Then, by 
definition of (v.a)g, all occurrences of a within g are overlined in (va)g (to 
denote the name extrusion). 

In rule (Com), the annotated tuples are “complementary” and can be 
merged, by merging both the link chains and the two tuples. Note that we 
leave implicit the side condition sges’g' # |, because the sges’g’ annotates 
the label of the conclusion transition. In addition, rule (Com) checks that 
(i) the extruded names of one process do not clash with the free names 
of the other process (like in ordinary z-calculus), and finally that (ii) the 
communication is not completed yet (se s’ not solid). 

Rule (Close) checks differ from the (Com) one, because (Close) is 
applicable only when the communication has been fully completed and 
cannot be further extended (se s’ is solid). In this case we must close, 
i.e. put back the restriction of all names extruded in the communication 
((v ex(geg’))). Still, we need to make sure that geg’ has no unresolved input, 
ie. that all requested values have been issued (ge g’ is ground). Moreover, 
the observed label is just s es’, as explained before. This is similar to the 
m-calculus mechanism, according to which the synchronisation of e.g. a(x) 
and a(x) yields 7 and not T(z). 

While rules (Lsum), (Rsum), (Rec) and (Struct) are straightforward, 
rules (Lpar) and (Rpar) need just to check that extruded names of one 
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process do not clash with free names of the other process (like in ordinary 
m-calculus). 

For some simple examples of the LTS semantics of link-calculus, we 
refer to [4]. Here, we just point out with an example the way in which the 
merge of link chains is performed. We recall that the definition of the e 
operator requires that the two involved link chains must be of the same 
length, and that, by rule (Act) we can arbitrarily add virtual links, in any 
position. Let s; =7\p, s2 =°\c, and s3 = °\, be three links, then there are 
three possible combinations for applying the e operator: 


e we can first define s, = 7\p*\s*\s, 8) = *\x?\c*\s, and compute 
84 = sesh = 7\,°\c*\e; then we can define s = *\,*\x°\,, and, 
finally, we can compute s = 4 s3 =7\,°\,°\,; or 


ee 


e we can first define si = 7\y*\.*\x, 53 = *\x*\«°\r, and compute 
84 = ses, =7\p*\«°\7; then we can define si, = *\,°\-*\x, and finally 
we can compute s = sh esq4=7\p,"\-°\,; or 


e we can first define 54 = *\,*\«°\7, 84 = *\s°\c*\«, and compute 
84 = ses, =*\,°\c°\7; then we can define s, = 7\y*\,*\«; and finally 
we can compute s = si es4=7\p?\¢°\z. 


We will see in the next sections many medium-sized examples of link- 
calculus at work for modelling biological mechanisms and for encoding Brane 
Calculi processes. 


4 Back to the Receptor-mediated Endocytosis Ex- 
ample 


Once the calculus has been formally introduced, we can revisit and complete 
our example, introduced in Section 2, based on the Receptor-mediated 
Endocytosis mechanism. 

To formally derive the first transition of the system, we start by deriv- 
ing the first part of the synchronisation between the LIGANDyop and the 
membrane Mer Note that we have extended the link prefix with virtual 
links in order to execute the merging of the two links. 


R R 
Piha ° ‘eo PN aeeishean oe We = T \iaemnen veeeapton : 
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*\imemRec\x\+_t oe i acepioe vein 
R 


top top 
—. > 
LIGAN Diop 0 mie » Mier 
T\ memRec 
t \rien Ree \receptor (top,MeMout,MEMin cell) t 
LIGAN Dtop|M.ch, > OMe 


where the tuples of LIGAN Diop and of M ee Lex LOD; igi nee yon) 


Cl 
and (top, MEMout; Lmemins Leell) Match, and, as a result of their combination, 


we obtain the tuple (top, memout, MEMin, cell). The value top can be assigned 
to the variable xo), while the value memo, can be assigned to the variable 
Lmem- Recall that our semantics is in early style. 

Now, we proceed our derivation with the synchronisation with the AP-2,¢1 
molecule: 


memRecy * 


LIGAN Dagp| Mion “SEREEE\ecertor\2\+ 


cell 


t” *\ *K\ K\ ap2 my 
> OL Mtr APO 55 \e\e\SP\7r t 0 


cell 


memRecy, * 


LIGAN Dtop| Mth [AP D.sh \memRes \recepter 


cell 


\ap2\ . (top,meMout,MEMin ,cell 


’, o|M'?? 10 


cell 


where the tuple of the first premise, (top, MeMout, MEMin, cell) matches with 
the tuple of the second premise, (top, MEMout, MEMin, Cell), and, as a result 
of their combination, we obtain the tuple (top, memout, MeMin, cell). 

The derivation of the transition can be concluded by the interaction 
with the RECEPTOR ce: 


Ne \eeeeP ee Na (top,meMout,MEMin cell 


RECEPTOReeii 5 LIG-REC-AP2 ei 


that leads us to the final transition, where the tuple of the first premise, 
(top, MEMout, MEMin, cell), matches with the tuple of the second premise 
(top, MEMout, MEMin, cell), and, as a result of their combination, we obtain 
the tuple (top, MeMout, MEMin, cell): 


LIGAN Dtop|M*2? |AP-2cett| RECEPTOReett 


cell 


{ Themen Re a a oN 
memRec \receptor \ap2 \T 


0|Mi2h |0|LIG-REC-AP2ce1i 

We can now pass to model the second part of the mechanism. We 
keep the same modelling style used for the first part, that allow us to 
model the creation of a vesicle membrane and its fusion with the endosome 
membrane. As a consequence, we focus our attention on the last two 
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MES pee X (ta ) 
+ 
o comple (cell, y) AX 
(the cell membrane) 


COMPLEX cen & ©°?!€?\ (a, cell) .(V,<2!!| COMPLEX ps1) 


vsl 


(within the cell) 


cell A 1 
Vosl =P \endoaanie (Lust; Lnew)-O 


(the vesicle membrane) 


COMPLEX ys1 & 7\ys1(vsl, Znew)-(AP-2cen |CLATHRIN cet |LIG-REC new) 


(within the vesicle) 


ENDOSOMEC! A CURIE sats new). NEWeell 


endo new 


(within the cell) 


Figure 4: Encoding of the Second Part of the Receptor-mediated Endocytosis 
Example. 
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biological interactions needed for the formation of the vesicle (V,°4!’) and for 
the releasing of the AP—2 and CLATHRIN molecules, as shown in Figure 5. 


In Figure 4, we give the link-calculus complete code of the example, 
where we follow the same modelling style used in Section 2. Here, we model 
the formation of the vesicle as the result of a further interaction between the 
COMPLEX,,1; (composed of the LIGAND, the RECEPTOR, and the two 
AP—2, and CLATHRIN proteins) and the cell membrane. Then, we model 
the uncoating of the vesicle and its fusion with the ENDOSOME as a single 
interaction between the membrane vesicle, the COMPLEX,,; (i.e. the same 
complex located inside the vesicle), and the ENDOSOME. This last reaction 
releases the AP—2 and the CLATHRIN proteins in the cell, and the fusion 
between the vesicle and the endosome gives rise to a new membrane inside 


the cell: NEW°e! 


new 


Clathrin-dependent endocytosis 


receptor formation of 
N clathrin-coated vesicles 


AP-2 
Cee cae <— 
4, clathrin \ 

uncoating 

AP-2 complex £2 @ 
a-adaptin 
62-adaptin «% al 
p12-chain fusion \ 


o62-chain 


Figure 5: The Receptor-mediated Endocytosis, Second Part. 
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M‘2? |COMPLEX,.u|ENDOSOMES),, 


cell 
complex 
1 TY omen 


M*°?| COMPLEX ysi| V2! |ENDOSOME“ 


cell vsl 


ee dea icbareactinag 
usl \endosome \T 


M°°? |LIG-RECnew|AP-2ceu| CLATHRIN cei |OINEW oe! 


cell 


Link-calculus offers a natural and simple encoding of the Receptor- 
mediated Endocytosis mechanism, thus encouraging us in investigating the 
possible use of link-calculus in directly modelling other biological phenom- 
ena. 


5 An Overview on Brane Calculi 


The Brane Calculi [8] are a family of calculi for describing biological mem- 
brane interactions. Membranes are not only containers, by constituting 
the boundaries of compartments, but they are also active entities: they 
represent the place where processes are located and interact with the sur- 
rounding environment. In this section, we briefly recall the Mate/Bud/Drip 
(MBD) version of the Brane Calculi [8], in order to show its encoding in the 
link-calculus, in the next section. 

As described by the following syntax, a membrane system consists of 
nested membranes, where each membrane is associated with a membrane 
process, that describes its interaction capabilities, in terms of possible 
membrane transformations. 


PQs= | Poe |!P|e«(P) systems II 
acr2= 0/7 |1c]ac membrane processes 43 
a,bu=  maten | mate} | budn | bud£(p) | drip(p) MBD actions Epp 


The basic structure of a system consists of (sub-)system composition, repre- 
sented by the monoidal operator o (associative, commutative and with o as 
neutral element). Replication ! represents the composition of an unbounded 
number of replicas of the same system or membrane process. The system 
o(P) is a membrane with content P and interaction capabilities represented 
by the process o. 
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(Par (Brane) (Struct) 
P+Q P+Q PEP APH3V A Q=Q 
PoRSQoR o(P) > a(Q) P-Q 


(Mate) maten.c|oo(P) ° matex.t|To(Q) > o|o0|T|70(P © Q) 


(Bud) bud} (p).t|T0(budn.o|ao(P) o Q) — plc|ao(P)) o T\70(Q) 


(Drip) drip(p).c|o0(P) > ple) ° aloo(P) 


Table 3: Reduction Semantics for Brane Calculi. 


The MBD set of actions are inspired by the possible spatial transforma- 
tions of biological membranes, such as membrane fusion and splitting. The 
former is modelled by the mating operation, the latter both by budding, that 
consists in splitting off exactly one internal membrane, and dripping, that 
consists in splitting off one empty membrane. 


Membrane processes o consist of the empty process 0, the parallel 
composition of two processes, represented by the monoidal operator | with 
0 as neutral element, the replication of a process and the prefixes of a 
process by an MBD action a. Actions for mating (mate,) and budding 
(bud,) have the corresponding co-actions (mate+ and bud-, respectively) 
to synchronise with. Here n, which identifies a pair of complementary 
actions, is taken by a countable set A of names. The actions bud‘ (p) and 
drip(p) are equipped with an argument p that represent a process, i.e. the 
actions associated with the membrane created when performing budding and 
dripping interactions, respectively. As usual, terms of the Brane Calculus 
can be rearranged according to a structural congruence relation =, whose 
rules (i.e. monoidality of operators o and |) are omitted for brevity (see [8]). 
The semantics of the calculus is given by the reduction rules in Table 3. 
Besides the standard reduction rule for congruence (Struct), and contextual 
propagation of reductions across parallel composition (Par) and membrane 
nesting (Brane) (cf. the upper part of Table 3), there are the axioms specific 
of MBD (cf. lower part of Table 3). 
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Rule (Mate) models the fusion of two parallel membranes. The effect 
of a mate synchronisation is twofold: (i) it creates a new membrane; and (ii) 
it dissolves the two existing ones. In the rule (Bud) a membrane expels a 
child membrane o|o9(P)) and surrounds it with a new membrane, associated 
with the process p, where p is the parameter of the co-action bud. Finally, 
in the rule (Drip), a membrane creates a new empty membrane associated 
with the process p, which is the parameter of dripn. 


6 From Membranes to Links 


6.1 The MBD Brane Calculus 


The encoding of the MBD Brane Calculus into the link-calculus follows the 
same style of the one provided for Mobile Ambients in [4]. The encoding 
relies on the idea that each interaction not only involves the two processes 
triggering the interaction, but it also has an impact on the membranes and 
processes of the context, that result reconfigured, after the reaction. To take 
the implicitly multiparty nature of interactions into account, the encoding 
offers a link-calculus counterpart for each actor involved. 

More precisely, the rule (Mate) requires a four-party interaction (at 
least), involving: 1) the membrane process with the capability mate,; 2) the 
membrane associated with this process; 3) the membrane process with the 
capability mate; 4) the membrane associated with this process. The rule 
has a double effect: (i) it creates a new membrane, and (ii) it dissolves the two 
existing ones. When the two membranes are dissolved, their contents must 
be relocated, which may consist of an unbounded (and not known a priori) 
number of parallel processes, all therefore involved in the interaction! As in 
[4], we adopt here a syntactic solution with no semantic impact on the rest: 
we replace the membranes with some sort of blind forwarders (reminiscent 
of those in [18]) that leave their contents unaware of the deletion. 

The rule (Bud) requires four-party interaction as well, while the rule 
(Drip) is simpler, as it only involves a membrane process and its correspond- 
ing membrane. 


Membranes with Brackets Forwarders impact on the encodings that 
have to deal with them, as they represent a marker (a side effect) of when 
and where a membrane has been dissolved. Since in the Brane code, no 
similar side effect is generated, the correspondence between MBD processes 
and their encodings is not direct. 
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To simplify our presentation, we directly introduce forwarders in the 
syntax of MBD, with no effect whatsoever on the semantics. We only need 
to slightly modify the syntax with the possibility to enclose a system P or a 
membrane process o within a pair of parentheses or brackets, as in [4]. 


Pease P| G2 ><|| [a] 


To make the presence of parentheses inessential with respect to the behaviour 
of the process, we introduce the additional structural congruence axioms: 


P=Q= [P| =([Q| g=7T => [ol =[7] 


Finally, we define the notion of passive context, both for systems (C), and for 
processes (S). Passive contexts allow us to adjust the basic reduction rules 
to deal with the presence of an arbitrary number of balanced parentheses in 
a monoidal composition. 


C,DiS=| |C) | Car | Pel S,T::=- | [S} | Slo | ofS 


and write C(P) to denote the system obtained by replacing the hole - in C 
with P. Similarly, we write S(o) to denote the process obtained by replacing 
the hole - in S with o. Thus we add the following suitable axiom: 


(Mate) C(S(maten.o)(P)) o D(T(mate;,.r)(Q)) + C([S(o)]|[T(7)]([P] ° [Q])) © 


(Bud) — T(bud;,(p).7)(C(S(budn.c)(P))) > pllS(o)| (PI) ¢ T(r)(C(o)) 


(Drip) S(dripn(p)-0)(P) > ple) oS(o)(P) 


Furthermore, we introduce the following reduction rule to propagate reduc- 


tions across contexts. 
P+Q 


[P] > [Q] 


In the above rules, we introduce only the contexts involved by the 
modifications due to the application of the rules. Note that, thanks to the 
use of contexts, we have eliminated some siblings of the interacting systems 
and membrane processes, which were necessary in Table 3. More precisely, 
oo and 79 of the original rules are now subsumed by S and T, respectively 
(in all rules), while the system Q of the original rule (Bud) is now subsumed 
by C. 


(Brac) 
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Structural Encoding The encoding function is defined in Figure 6, and 
works similarly to the ones in [7, 4]: 


e The hierarchy of membranes is rendered as a flat system of parallel 
link-calculus processes, where the inclusion relation is recorded by 
interaction channels. 


e The encoding | P |] of a process P is parametric with respect to a given 
family of channels m that represents the enclosing membrane or, the 
topmost compartment top in which the membranes reside. Informally, 
each process resides in some “location”, to which it addresses all its 
requests about the next actions to perform. 


e Each membrane o(P) is encoded by a link-calculus process of the 
form M? ;, in parallel with the process | a]; that encodes the associated 
membrane processes o, and with the one [| P],;, encoding the included 
process P. Intuitively, M? ; represents the “control or management unit” 
of the membrane, a sort of interface that offers a family of channels, 
or ports, for each kind of action capabilities. Families of channels are 
indeed indexed by subscripts mate, matet, bud, bud+ and drip. In 
the following, such families are called locations (in Figure 6, families 
are written for brevity as m,m’,... (and the like)). The links of these 
channels are composable with the corresponding channels of: (i) the 
encoding of the actual membrane process o, that represent the active 
capabilities of the process, (ii) the encoding of the parent membrane, 
and (iii) the encoding of the included process P. More precisely, in 

MP jr tis the family of channels shared with | c]j;, j is the family of 

channels shared with [| P];, while p is the family of channels shared 

with the parent membrane. 


e After a fusion, the membrane hierarchy changes: the sub-membranes 
of the two fused membranes become siblings, and, as a consequence, 
their interaction capabilities change. On the encoding side, this implies 
that locations must be updated and that the interactions relative to 
the processes originally included in the dissolved membrane must be 
redirected. To perform this task, the translated processes make use of 
the forwarding mechanism every time a membrane is dissolved. The 
forwarding activity is managed by processes of the form Fwd?, , located 
between the dissolved membranes and the parent ones. 
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[o]m = 0 [P°Q]m = [P ]m|1Q Jim 
[!P]m = recX.([P]mIX) [[Pl]m = (ym’)(Fwdy, | P]m*) 
[oUP) Im = (Yumi, me2)(MR mg | Lo Im | P Ime) 

[0] m = 0 [o|T]m = [o ml 7 ]m 
[!o]m = recX.([o]mlX)  [lol]m = (m')(Fwdty | Lo]m) 


[ maten.o im 

[ mates .o]m 
[ budy-o Jin 
[bud (p).o]m 
[arip(o).o lm 


7 \mmate (Ms 71) £2).[ o]m 
™Mmatet Veta £2). 0 m 
e eaied (n, T1, x2). 0 Im 
™ouat \-(n, 01, 02)-([p Jes | ]m) 


*\marip (€1)-(LP Des | 7 Im) 


[> [I> [I> [I> I> 


MP. £* recx. (vm1,m2)'™*** \p.a4- (AME, mM1,M2).( Fwd)" Fwd}? | X[m1/i,m2/J]) 


a,j 
+ Pmate\ name, £1, £2).(Fwd;" | Fwd;*) 


matet ( 


af toud\ ug (NAME, 1, £2).(X[x2/pI) 
(vmy,mo)"4\;,_, (name, m1, m2).(X [my /i,m2/J] |X) 


+ (vm1,m2)'4"?\,(m1).(X [mi /i,m2/j] | X) 


Fwd? = recX.( ™mete\p (mame, £1,%2).X + Pmate+\m — , (name, x1, %2).X 


atet 
a de re (name, £1, ©2).X cee \my ua (name, £1, a2).X 


aes eat (x1) X 


Figure 6: Structural encoding of Brane in link-calculus. 


We only comment on the translation of the reductions (Mate), (Bud) 
and (Drip). 

As said before, the mate synchronisation is a four-party one. For the 
sake of clarity, we focus the discussion on the Brane code below, illustrated 
in Figure 7, where we tag the involved membranes and brackets by a, b,c 
(see below) to refer them in the text. 


maten.7|o(P)a o mate; -T|70(Q)» > [oloolal[t|Tolol[Pla © [Qla)e 
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membrane 


Figure 7: Illustration of the Mate Synchronisation Encoding: the Source 
Process. 


First, we note that [| maten.c|co(P)a o mate+.|79(Q)s Jim amounts to 


(v a1, a2, bi, b2)( Mat aglll maten-o Jail all P das | 


Myr poll mater 7 Jo, IL 7 Jos IQ Too ) 


e The mate,.c party: the process | mate,.c ja, provides the link prefix 
\ imate (Ms £1, 2) that requires the synchronisation with the corre- 
sponding link, offered by the translation M7’) ,, of its membrane. Here 
the message n is sent, and the names x1, X2 are received (but they are 
irrelevant). 


e The party M;’' ,, for the membrane a: the translation of the membrane 


a holding the mate, operation is: 


recX. (17:61, 62)" \ rmare (NAME, C1, Co).(Fwde | Fwd? | X[e;/ai])+... 
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Figure 8: Illustration of the Mate Synchronisation Encoding: the Target 
Process. 


where i € [1,2] and c1,c2 are two freshly created locations to be 
communicated to all the other parties, X[c,/a1,c2/az| can be just 
read as Mf’... (ie. it represents the new membrane c), the link 
@1mate\, .. is needed for the synchronisation with the mate,.o party, 
and two new forwarders, Fwd(\ | Fwd? are created for redirecting 
the communications with components that are still located at the 
updated locations a, and az (because such components do not actively 
participate to the interaction). Note that the variable name is used to 
receive the parameter n of the (Mate) reduction, even if it is not used 


in the continuation. 


The party M;",. for the membrane b: the translation of the membrane 
b holding the mate+ operation is even simpler: 


PeGX 2 so a nately (name, £1, ©2).(Fwdy" | Fud;;) +... 


Imatet 


A Flat Process Calculus for Nested Membrane Interactions 117 


where the link "™ete\,, is for the synchronisation with the trans- 


lation of mate+.r party, and two new forwarders F’ wd, | F wd; are 
created to redirect the synchronisations on the old locations b; and b2 
to the received locations x; and x2 (they will correspond to the freshly 
created locations c,, cz extruded by the party M’” ,.). 


a1,a2 


matet 


e The mate;.r party: the translation of mate+.r is similar to that of 
maten.o. It provides the prefix Pimate\,(n,x1, 22), where Dymatet is 
for the synchronisation with the encoding of its membrane, the message 
n is sent for the matching with the one sent by the mate,.o party, and 
the received names 71, 29 are still irrelevant. 


Putting all pieces together, we can observe that the four-party synchro- 
nisation forms the following link chain, once reached the agreement on the 
associated tuple (n, G, €) 


ees Mmate\ imate! \ 
Almate \Mmate b, T 


mate 


where the freshly created names c),co are extruded to the other parties. 
We can then observe the correspondence of the resulting process with the 
encoding of the target system [o|oo]a|[T|To]o({Pla © [Q]o)c, illustrated in 
Figure 8, by noting that | [o|oo]a|[T|To]o([Pla © [Q]b)c |m amounts to 


(v a1, a2, b1, be, c1,¢2)( Moye, | Loloolaill P Ja2| Pwd) | Pwde, 


[7170 ]o IL Too! Fwd, | Fwd, ) 


As, in general, we would need to consider the presence of contexts C 
and S (i.e. of brackets in the syntax), the four-party synchronisation can 
become a larger multiparty synchronisation, where an unbounded number 
of forwarders is involved. The same holds for the cases we discuss next. 

Also, the translation of the bud, operation involves four entities. In 
this case no forwarder will be introduced. The Brane code we refer to is: 


bud;,(p).7|To(budn.o|o0(P)a° Qs + plo|oolP)abc © T!TlQ)e 


where we still distinguish the two involved membranes a and 8, holding the 
two co-actions bud, and bud+(p), respectively. 
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The bud,.o party is translated as 7\q,,,,4(N, 1,2)... ]a,- The link 
*\ aimma 18 for the synchronisation with the translation M® o of its 
membrane, the parameter n plays the same role as before, and the 
received locations x1, £2 are irrelevant. 


The party MP. for representing the membrane a is: 


recX. +++ “lbud\, | (name, £1, £2).(X[x2/be]) +... 


where X[x2/b2| can be just read as M7? (ie. it represents the mem- 
brane a after the relocation in the new membrane c), and the link 
“1bud\y..,q 18 needed for the synchronisation with the budy.o party. 
The variable name is used to receive the parameter n of the (Bud) 
reduction, even if it is not used in the continuation. The variable x, is 


irrelevant. 


The party M;",, for representing the membrane } is: 


recX. ---+(v cy, oa) May: name, C1, C2).(X[e1/b1, co/b2] | X)+... 


budt ( 


where two fresh locations c,,c2 and a new membrane X [c;/b1, c2/b2] 
(to be read as M2” ...) are created to represent the new membrane c, 
and the link PAN 
of bud (p).r party. Note that locations c;,c2 are extruded in the com- 
munication, as they will be needed by the other parties to (re)position 


p and the membrane a. 


, is for the synchronisation with the translation 


The bud (p).r party is translated as Piua+\,(n, 21, ¢2).([ per |[ Jo, )- 
The link iowa \, is for the synchronisation with the translation of its 
membrane M, ane the parameter n plays the usual role, and the received 
location x1 is used to locate p (only the location x2 is irrelevant). Note 
that x, will be substituted by the freshly created location cy. 


The four-party synchronisation forms the link chain below, once reached 


the agreement on the associated tuple (n, , @) 


7\ Qbud bobud\ Proud \ 
@1bud \bobud \Oyy,,q-1 \7 


where the freshly created names c;, cg are extruded to the other parties. 
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Finally, we can discuss the translation of the action drip, that, since it 
involves only one membrane, is the simplest one: 


drip(p).o|70(P)a > plo)e ° o|a0(P)a- 


The drip(p).o party is translated as 7\a,,,;,(€1)-([ ]ai|I 7 Ja,), while the 


party M;'a. for representing the membrane a is: 


recX. +--+ (Vv b1, bg) 4rp\ (by) .(X [by /a1, bz /a2] | X) 


where the fresh locations b;, bz and the process X [b;/a1, b2/a2] (to be read 
as My ,) represent the fresh membrane b to be used for locating p. The 
two-party synchronisation forms the link chain below, once reached the 
agreement on the associated tuple (b;) 


ae ‘ 
G1drip \T 


Example We will now show a fully detailed example that illustrates the 
encoding of a MBD process in link-calculus, and also two of its possible 
derivations. Let P = mate,.drip(p) | o(r(R)) ° matet.6(Q) be a process, 
where two membranes merge, and then a drip reaction creates a new empty 
membrane. 


P > P, = drip(p)|o|6(7(R) o Q) (mate) 
Pi P2 = pl) oold(r(R) oQ) (drip) 


In the following encoding, we colour the parts involved in the mate transition, 
and we call, for brevity, 01, 72, 63 the translations of the corresponding brane 
actions, where the index is the number of the membranes that initially 
include them. The encoding is | P Jiop = Prinx|Pling, Where 


Prine = (Y'71, 81)(Min? 5, |” \mimare (% @1)£2)-7\mrarep Y1)-(LP ]y,!o1) (membrane m1) 


(v m2, 82)(Man 50 |72|[ R]s.)) (membrane mz in mj) 
Pling = (¥ M3, 83)(Mig? gg | ™smatet \-(n, £1, £2)-03 | [| Q]s5) (membrane 73) 


For simplicity, we just give a sketch of the derivation of the two operations. 


t 7. = % 
T\mimate\,..° | \t\x(nm! *\z\ formate Mamatet \ , (n m/s! 


3’) i) 
; r 3mate ! 
> Qrink Phink > Qiink 


T\T\ foPmate\ 7 \ 


1 °Pmate i 
Prinn| Phin Prasite 


Prink 
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In Phy, = (vm',s’)(Qrink|Q ink), we colour the parts involved in the drip 
transition, where also the forwarder process is involved, as the membrane m, is no 
more active: 


Phe = [Pilon = (Ym! 8") (vm, 81)(Fude, | Pwedl, | MY? 4 |\ mrarsy (Us) Planter) | 
(v m3, 82)(M&, ay 1721 [Rl sa)) 
(v ms, 83)(Fwd?", | Fwd, |63 | [@]ss)) 


*\ 1drip * u” , 

\ * \in! . \x(m™) x\x\ "drip \ (m!’) 
pla drip gia pib * «Ar qle 
link link link — Os @ Dink 


7 
™1 drip, "dri ” 
RNP \ ah!) T\ mz gigy Ve Ve Gm!) 
: Lee 
la 1b drip la 1b le ldrip le 
Pritxl Pink > Qik! Qiink Pink Qik 


f 
M1 drip, drip 
2 
1 1b )pl \madrip mney 1 1b iol 
a c a e 
Pring | Pijox | Pring Q1 ink! @Qiink | Qiink 


mi 
PRN pe oe Na 


1 drip 2. 
Prink Prinx 


where P?... = [ P2 top is equivalent to 


(um, s') ((um"',s”)((v m1, 1)(FPudm, | Fwds, | Mie? | Ma? go | Lem) 101) | 
(v ma, 82)(Mr, s. 1721 R]s2)) | 
(v ms, 83)(Fudy, | Pwds, | d3|[Q]s.)) 

The result of the above transitions is that membranes m1 and m3 have 
been substituted by the new merged membrane m’ lying on the top position, 
and the four processes F’ wdm , F wds, vin wd .F ws, keep track of this event. 
Another new membrane m”, lying on the top position, is introduced as the 
result of the drip operation. 


Operational Correspondence We can state an operational correspon- 
dence result to show the correctness of the proposed encoding. 

Intuitively, the idea is to define an encoding assigning to any MBD 
system P a corresponding link-calculus process [ P ]. 


e for any reduction P > P’ there is a step | P] > | P’]; 


e and vice versa, for any silent step | P] — Q there is a MBD system 
P’ such that Q = [P’] and P > P’. 
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Unfortunately, a direct encoding has to deal with the presence of forwarders, 
so that in general: 


e for any reduction P > P’ we can find a step [P] — Q but Q can 
differ from |] P’] because of the presence of forwarders; 


e and vice versa, for any silent step | P ] — Q we can find a MBD system 
P’ such that P + P’ but, again, Q can differ from [ P’] because of 
the presence of forwarders. 


We obtain a correspondence based on the forwarders introduced in the syntax 
of MBD, with no effect whatsoever on the semantics and expressiveness, and 
that allows us to recover the stronger correspondence result sketched above, 
with exact matching between the MBD reductions and the silent steps of 
the link-calculus (modulo some standard structural laws imposed by the 
MBD structural congruence). 


Proposition 1 The encoding in Fig. 6 is well-defined, in the sense that if 
P=@Q then[P] =[Q]. 


We can filter out non solid transitions (representing partial interactions) 
of encoded processes by letting: [| P] = (Vv tmate; timate» tbuds thud» tarip) | P lt 

We say that a link chain s is silent if it consists of 7 actions only, and 
write P > Q if P > Q for some silent s. (Note that any silent link chain is 
solid.) 


Lemma 1 /f[P] > Q, then s is silent. 


Sketch of the proof: The proof is by cases on the type of capability 
simulated by the encoded process [| P] (mate, bud, or drip,). The proof 
of each case is, in turn, by induction on the length of the derivation of the 
transition. 


Theorem 1 Let P be a MBD system, then P — P' if and only if there 
exists Q such that |P] > Q, andQ=[P’]. 


Sketch of the proof: The proof for the only if part is by induction on the 
length of the derivation of the reduction P — P’; the if part can be proved 
by cases on the type of the capability simulated by the encoded process. 
Here, some explanations about the derivation of a silent move: [| P] > Q 
are in order. First, the silent move s must have been obtained as a chain 
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of link prefixes that come equipped with compatible tuples, and that are 
combined along restricted channels. Of course, at the left and right hand 
sides of this link chain, we need two link prefixes presenting a 7 on its left 
and on its right, respectively. Then, thanks to rule (Close), in Table 2, we 
can derive a link chain s where only 7 appears and the tuple of values has 
been discarded. 

Now, we detail the construction of the link chain. At the left hand side 
of the link chain, a link offering a silent action 7 is required. By the rules 
in Figure 6, the links offering an action 7 on their left are: 7\a,,a:¢5 ”\apua> 
and 7\q,,;, deriving from the encoding of the prefixes maten, budn, dripn. 
At the right hand side of the link chain, a silent action 7 is also required. 
By the rules in Figure 6, the only links offering a 7 action on their right are: 
bmatet\., buat \, and %rip\_, The first two links derive from the encoding 
of the prefixes mate‘, and budz; the last one, %"»\,, derives from the 
encoding of the membrane holding the prefix drip,. Back to the encoding 
of the capability mate, two membranes are needed to pass from the channel 
Gmate to the channel Mmate, and from b,,44¢1 tO Mmate: on their surface, the 
two membranes host the mate, and the mate; prefixes, respectively. As 
fas as the encoding of the capability bud is concerned, two membranes are 
needed to pass from the channel apy,g to the channel mpyq, and from Mpxyq 
to dé 4: on their surface, the two membranes host the bud, and the bud- 
prefixes, respectively. The case of the capability drip is simpler, because it 
does not require additional links. Forwarders can provide links to the chains, 
but without changing the type of the channels. Participants to the action 
can be arranged in parallel, according to their position in s, by using the 
link-calculus structural congruence rules. By induction on the proof of the 
derivation, we can rebuild the tree of the bracketed membranes involved in 
the reduction. 


6.2 The PEP Brane Calculus 


We now extend our framework, in order to include also the Phago/ Exo/Pino 
(PEP) version of the Brane Calculus, where the three new reactions represent 
the operations that model endocytosis and exocytosis. The first indicates 
the process of incorporating external material into a cell, by engulfing it 
with the cell membrane, while the second one indicates the reverse process. 
Endocytosys is rendered by two more basic operations: phagocytosis (phago), 
that consists in engulfing just one external membrane, and pinocytosis (pino), 
consists in engulfing zero external membranes. Exocytosis is instead denoted 
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a,bz= —... | phagon | phago+(p) | exon | exox | pino(p) 


(Phago) phagon.o|oo(P) ° phago;,(p).T|70(Q) > T|Toleo|oo(P)) o Q) 
(Exo) exo .t|T9exon.o|a0(P) o Q) > Po aloo|T|70(Q) 
(Pino)  pino.(p).c|ao(P) > oloo(p(e) o P) 


[ phagon.o im 

[ phago;, (p)-0 |m 
[ pino(p).o ]m 

[ ex0n-0 Im 
[exot.o ]m 


*\ mphg (Mr £1) £2). 7 Im 

prot \,-(n, £1, £)-([e]e, | Lo ]m) 
\mpino £1) (Lp Jar | Lo ]m) 

*\ meno N, £1, £2)-[ 7 Im 

enol \e Wt Bi, £).[o]m 


I> [I> [I> [> I> 


M?, 2 recX.s +o + ‘pha \y 4g (RaMe, £1, ©2).X [x2/p| 
+ (vm, M2)PPhs roee (name, m1,m2).( Mgr ms |X) 
+ ‘ero\,.,, (name, x1, t).(Fwd;* | Fwd}) 
+ Jero\; |, (name, i, p).X 
+ (ym1, M2)*Pine\, (m1). (M™ _,,|X) 


™m1,™2 


Table 4: The Phago/Exo/Pino (PEP) Syntax, Semantics and Encoding. 


by exo. The extended syntax and semantics are in the upper part of Table 4, 
while the corresponding encoding is in the lower part of the table. 


As usual, the membranes are encoded as separate entities, and each 
synchronisation involves the membrane and one or two prefixes, correspond- 
ing to the operators pino(p), or phagon and phagos(p), or exon and exo. 
For the operators pino and phago, a new membrane is created and extrusion 
is exploited to relocate pieces of code: actually some biological elements are 
forced to change their location, as they are sent to the new membrane. We 
omit the obvious extension of the definition of forwarders and the detailed 
discussion of the various cases, as the technical development is analogous to 
the MBD case. 
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6.3. Comparing the MBD and PEP encodings 


As discussed by Cardelli, the MBD set of membrane capabilities (mate, bud, 
drip) can be derived from the PEP membrane primitives (phago, exo, pino). 
Nevertheless, since they model particular biological phenomena, it is worth 
having them as primitives as well. 

For instance, the mate interaction can be obtained by a suitable combi- 
nation of phago and exo, as illustrated below. 


A 
mater = phagon.exon 


mate+ * phago+(exo+.exo,).exo>y 
n n n n n 


It would be interesting to investigate the connection between the MBD 
encoding of membrane primitives and the similar PEP one into the link- 
calculus. To this aim, we try to compare the encoding of the MBD process 


Pusp maten.o|ao(P) o mate+.r|79(Q) 

[Pusp]t = (vnl,n2,m1,m2) 
Mrinal” \ntmate (£1, £2).[ 0 Ini|[ 70 Jnill P Incl 
Mint mal *ete\7(r, 21, £2).[ 7 alll 70 Jal @ Ime 


and the encoding of its corresponding PEP version, obtained by using 
the above combination of phago and exo. 


Pprp 


phagon.€£0n0\o0(P) ° phagox (exoy exp") 020-3 .T|T0(Q) 
[Peep le 


(v nl, n2,m1,m2) 


T\ nlpng (£1, £2).7\ nto. (n’, £1, £2). 0 Intl 0 Intl Mir mall P Ina! 


mora\-(n, 21, 22).(Lexod exon” Jor |2\,(n", v1, £2).[ 7 ]mi)| 
[ 70 Imi |Mint mall Q Ime 


We show the transition of | Pyygp ] corresponding to the mate interac- 
tion, and then the three transitions of | Ppgp ] corresponding to the three 
interactions needed for simulating the mate (one phago, and two exo). Fi- 
nally we compare the two results. To help the intuition, we leave the channel 
names in the link chains and the tuples of values in the transition labels, 
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although the semantics of the language would remove them. 


[Pusplt = (vni,n2,m1,m2) 
Mya nl” \nimate (Ms £1, £2). 0 nil 70 Jnill P Jnal 
aL 
Myvi mmol" mate\- (1, £1, £2).[ 7 Jill 70 Jill Q m2 


mate L T\ Rl mate tmate\™ dinate \ dn ae 72) 


N1mate \tmate \M1mate 


(vnl,n2,m1,m2, newl1, new2) 
[oJnill 0 Ini|Fwdner| Pudge Mi ewt,newalll P Ina 
[7 Jill 70 Imaal Fwd? | Fwdis’?|[ Q Im2 


| Peep i = (vn1,n2, m1, m2) 
T\ nlpng (% £1, £2).7\nteo(n’, £1, 22).[ 0 Intl 0 nal My alll P Jnal 
ak 
ina\ (m1, 22) ([ eof exo Jeil"™\(" a 22).[7 Jt) 


[ 70 Ire] Mint mall Q Ime 


1 
phago dig Os nal tile te mero (n, new1, new2) 


(v nl, n2,m1,m2,new1, new2) 
\ nero (M's #1, £2). Lo Intl. Ina Merl P nal 
siete se) v1, £2). "\newhesa(M x1, £2)|[ 7 ]m1)| 


[ 70 Jil M, Lb cea VE Fatmall Q Jm2 
it 
exon’ | OMe eeu eeu (n’,new1, m2) 


(ynl1,n2, m1, m2, new, neue) 

[o Juall a0 Joi Pwdlst Pwd IEP Jnal 
\newhese (M xi, #2) 7 Jal 

[70 Jil M, aa aenae mi,m2ll Q Jon2 


” 1 2 1 
exon | \ newt e \n ae Veer, mI, t) 


(vnl,n2, m1, m2, new, view), 
[oIoalf oo Juul Fende gt Podrs ALP Jnalll 7 mall 70 Imi! 
Fwdy. |Fwd! owolM, mall Tine 


We can now compare the two final configurations (where the systems 
are rearranged in order to facilitate the comparison). 


aa 
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(MBD) 
(v nl, n2,m1,m2,new1, new2) 
[o Jnall oo Jni|Fwdng!|FPwdng?| Mt owt mewall P Inal 


nil 
[7 Jmtll 70 mi | Pwd |Fwdt gs’? |[ Q ]me 


(PEP) 
(v nl, n2,m1,m2,new1, new2) 
[oJnill 0 Ini Pwang | Pwdr? |Min mall P Inal 


ni 


[7 Jmll 70 mi | Pwd, |Fwd! wall Q me 


As it can be observed, in both cases we obtain the same hierarchical 

membrane structure. More precisely, in both cases we obtain a unique 
membrane, called Mjewinew2 in the MBD case, and Mf, 2 in the PEP 
one. The two membranes hold the same subprocesses and expose the same 
capabilities on their surfaces. Thus, we can say that we obtain the model of 
the same biological system. 
From the syntactical point of view, and looking at the different forwarder 
processes, we can instead observe that there is a different use of the restricted 
names in the two cases. This phenomenon is due to our encoding technique, 
which uses the forwarder processes to leave track of the dissolved membranes. 
As a consequence, the different membrane evolution histories in the two 
computations are recorded in different forwarder processes. This seems 
reasonable, from a biological point of view, because it reflects the different 
traces, produced in the two cases. Technically, this amounts to a slightly 
different parsing of the same term, due to a different insertion of brackets. 


6.4. Molecules and Molecular actions 


In Brane calculi [8], Cardelli also models small molecules that can easily cross 
or be transported across membranes, and whose movements are completely 
mediated by membranes. The new set of actions models the binding and 
the releasing of the molecules on both sides of membrane surfaces. 

The molecular extensions of the syntax and of the semantics of the Brane 
calculi is in Table 5. We slightly modify the original syntax by considering 
only the set M of molecules (ranged over by t, 1), i.e. we exclude molecular 
complexes. 

The encoding of the molecules and of the molecular reaction are in 
Figure 9. 


A Flat Process Calculus for Nested Membrane Interactions 127 


Systems P,Q s:= ...|t systems with molecules t,r € M 
Actions a,b s:=  ...|ti(te) 3 t3(t4) bind & release of molecules 


B&R t1, 0° ty (t2) = t3(t4).o|o0 (to ° P) — t30 a|oo(ta ° P) 


Table 5: Molecules and the Molecular Reaction. 


We slightly modify the encoding function style, to keep track of both 
names associated to a membrane, i.e. [ PJ], so that we initially have 
[P] = (vt) P]f. The only encoding rule that results modified is the one 
for membranes 


[odP) Im = mi, m2)(Min mo Lo Vm? [LP Tin3), 


while the other rules will keep the second membrane name as superscript for 
the encoding function [...]"”’, without any change. 


For the encoding of the B&R operation we distinguish four cases de- 
pending on the number and the position of the molecules that will disappear. 
We also make use of three names: in, out and inout (in bold in Figure 9 to 
point out that are special names) to record the fact that the molecule which 
is inside, or outside the membrane, or both the molecules, will disappear 
after the B&R operation, respectively. The two prefixes offered by the mem- 
brane M7, m, are equipped with the link \,,, that will connect the outside 
molecule with one of the previous prefixes encoding the B&R operation. 
Finally, a molecule t is rendered as a process Mol;” located in membrane 
m that can recursively fire two types of prefixes: one at the beginning of 
the link communication, \,,, and the other at the end ™\,; for each type 
of prefix, according to the encoding of the B&R operation, there are four 
possibilities, depending on the number and the position of the molecules 
that will disappear. 


As an example, we show the encoding of a molecular reaction, in which 
both the molecules disappear: p; © pi(p2) = o(¢)(p2 0 P|). Our encoding 
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(Vv 


[tI 


[ t1 (ta) =s t3(ta).o |? 


Me 


m™m1,m2 


Mol” 


(I> 


(> 


(> 


I> 


Mol;? 
i \ms (tr, ta, ts, t)-Lo Wie 
m1 \ mg (in, t1, t2, empty, ta).[.o 7"? 
( 
(i 


IN aa out, ty, to, ts, empty). o He 


ee Nine inout, th, ta, empty, empty). [ o ii 


recX. ... + 
TON (t1, ta, ts, ta). Min, m2 
+ 


™\ mi (type, ti, t2,t3,ta)-Mir ma 


recX.( 7\m(t, to, tz, t4).X [t/t] 


+ 7\m(in, t, to, empty, ta) X [ta/t] 


alt sees 


ut, t, t2, t3, empty) .0 


q 


m 
m 


rin, t1,t, empty, ta).0 


m 


t 
(i 
(o 
(i 
(ti, t, tg, ta).X [t3/t] 
(i 
( 
(i 


\ 
\ 
vn 
\ 


+ 7\m(inout, t, t2,empty, empty) .0 


7(out iti, t, ts, empty). X[t3/t| 
(inout, t;,t, empty, empty).0 


iftzs FoAty 
iftz =o 
ifty =o 


iftgz =o=ty 


Figure 9: The Encoding of the Molecules and of the Molecular Reaction. 


generates the following link-calculus code: 


t)(Mol;,, 


|(vmi,m2)( 


Mina |” \mz (inout, pi, p2, empty, empty) | Mol?) 


L 7\f\2\r2\, (inout, pi, p2, empty, empty) 


(vt, m2)(O |(v m1)(Min, m, |9| 0)) 


m1 ,™2 


where the transition label has been kept for clarity, and each link in the 
label has been offered by the subprocesses in the corresponding position, 
i.e. the first link, 7\;, has been offered by the first subprocess Mol! 


on. 


and so 


This further encoding shows how the link-calculus can easily simu- 


A Flat Process Calculus for Nested Membrane Interactions 129 


late different biological entities, considering also their locations, without 
introducing any new operator. 


7 Our encoding at work 


To show how our encoding works, we consider the example of the Semliki 
Forest virus, presented in [8]. Briefly: (1) a virus enters a cell via phagocy- 
tosis, then an endosome compartment merges with the virus, and the virus 
deposits its nucap particle (NUCAP in the code) into the cytoplasm, i.e. the 
liquid inside the cell, via an exocytosis (Figure 7, in [8]); (2) a series of 
molecular reactions reproducing the RNA virus replication follow (Section 
4.6, in [8]); (3) the virus reproduction cycle ends with two further steps 
(Figure 8, in [8], here reported as Figure 11). 


For the sake of simplicity, we only show the execution of the encoding of 
the part (1), which consists of membrane interactions, and the first reaction 
of part (2), which amounts to a molecular interaction. A nucap particle is 
defined as a membrane containing vRNA. The nucap surface can disassemble 
the nucap membrane, by pushing the molecule vRNA outside, in response 
to some trigger molecule found in the cytosol, (molecule dis-trg in the code). 


The Brane code is as follows: 


phagon.ex0,(NUCAP) o (virus) 
phagox (mate, ) (mates |exoz (o) o dis-trg) (cell) 
where 


NUCAP = dis _trg(vRNA) = vRNA(o) (vRNA) 


As there are no molecular interactions, we use the standard version of 
the encoding function that only uses one subscript name, i.e. | P],. The 
corresponding link-calculus code is as follow: 


t * 
Phagoy,s1-EX0ursi|Mjs1,urs2| Nucap| (virus) 
+L t 
Phagozeu|M eeu cett2| (cell membrane) 


ali alt cell2 cell2 A ‘ 
MatejorelEXOpgre1|Moaret,bare2|MOlgis ing (inside the cell) 
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Phagoyrs1 = "\vrslpng (n, #1, £2) 

Phagoz,, + cell gt \e(; Bis 29). matey |x: 

Nucap = nucapl\ » icap2(M1, Mg, dis_trg, vRNA, vRNA, empty) | 
MT heusape | Molen 

Mateparet 5 are nate! \r (n, #1, 2) 

EX0yrs1 = ON iets (n, v1, £2) 

ExOpre1 = Pareles \r(n, 71, £2) 


To help the intuition, we leave the channel names in the link chains in the 
transition labels, although the semantics of the language would remove them. 

In Figure 10, we show the execution of the link-calculus process. Ini- 
tially, there are a virus, whose membrane is encoded as VE si Aedes and a 
cell, M171 ce. Both lie next to each other, in the t¢ location. With the 
phago reaction, the virus enters the cell: the new membrane vesicle is cre- 


os ae cell2 : . 
ated inside the cell (Mel aicwieien)s and the virus membrane changes its 


location (Mea Then, the bare membrane, which models the endosome 
compartment inside the cell, merges with the vesicle. The effects of this 
operation are that (i) a new membrane new is created, and that (ii) the two 
membranes, vesicle and bare, are dissolved. As a side effect, four forwarders 


newl new2 newl new2 : 
Fudyesicter: HW yesicleg, LHWapare, and Fwdy\c5 are created to redirect any 


request from the old locations to the new ones. The last membrane reaction 
is the exo that allows the virus to exit the membrane new, and to dissolve 
itself, leaving proceses Fud’©!| Fud%'2. Please note that the forwarder 
Fwd"©" redirects the actions, previously located on the surface of the virus 
membrane, on the surface of the new membrane; this is a direct consequence 
of the behaviour of the action exocytosis, see Table 4 for the details of the 
Exo transition. The nucleocapsid VS i aaicne can interact with a molecule 
that lies in the cytoplasm (7 Ole. q) which trigger the disassembly of the 
nucleocapsid. Thus the mol interaction can take place. The result is that 
the molecule vRNA, initially inside the nuclecapsid membrane M Mas rev is 
now free in the cell Mol¢¢2 ,, and that the M Oleg has been “consumed”. 

Our encoding renders each biological location, i.e. each membrane, as a 
process that takes part in the multiparty interactions. Forwarder processes 
can record the history of the computation, giving the exact position of the 
membranes when they have been dissolved. This could be useful for an exact 
analysis of the dynamics of the computation. Alternatively, we can instead 


garbage these pieces of code, as in [7], by applying a renaming operator. 
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t . 

Phagoyrs1-EX0ure1| Moret vrsa|Nucap (virus) 
t 
Phagoceitt|Mccitt cettal (cell membrane) 
cell2 cell2 . . 
Mateparei bare2lEXOpare|Mecret,bare2| Molais trg (inside the cell) 
hy PART eka tore cen het \ 
phago vrslpng \tphg celll | ot f. 
t 

Meetti,ce2| (cell membrane) 
(mesiclel vestcle2\( \westdiel sj tia ees Wren setae (vesicle inside the cell) 
ExOvre1|MyrstsorsalNucap)| (virus inside the vesicle) 

cell2 cell2 . . 
Mateparer|EXOpare1|Mocret,bare2| Molais trg (inside the cell) 


‘ 1 
mate |7 vesiclel mate cell2mate\ PAPE lmate\ 


i. abe 
vesiclelmate \cell2mate \barelz |,. 


Moet cetg|(v vesiclel, vesicle2, new1, new2)( (cell membrane) 

Pudi ieer Fwd erica | Fwd pene | wd pines (dissolved membranes) 
Ex0yrs1|Mpesicle2 .|Nucap| (virus inside the new membrane) 
Exo pce Me unio, (new membrane and ) 


(dis-trg molecule ) 


exo |7 vrslero\ vesicle2exo new2era\NeWlero barel ene 
urslero \vesicle2exo \neW2ero \new1z,, \bareld,,, \ST 
Mt cetlz|(v vesiclel, vesicle2, new1, new2)( (cell membrane) 
mew 1 cell2 mew 1 new2 mew mew2 . 
Fwd | Pudee | Pudi | wd o|Fwdpes |Pudpes (dissolved membranes) 
meee \ncap? cOUE ta ieee NN eocaabey) |Miree pl.nucap2| (nuclecapsid membrane) 
cell2 cell2 
Mrewl,new2|M ol gis trg) (new membrane and) 


(dis-trig molecule) 


mol | T\ cell2\ urs2 sea en Sanka ieek 
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Figure 10: Execution of the Encoding of the Brane Phago, Mate, Exo and 
of a Molecular Interaction. 
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Figure 11: Viral Infection (highlighted part) and Reproduction. [Adapted 
from [8]. 


8 Conclusion 


We presented the application of the link-calculus to Systems Biology, show- 
ing that it is suitable for representing biological interactions. In this field, 
interactions vary on the number of elements, on their locations, and on their 
shapes. Here, we have shown that the link-calculus can easily take into 
account interactions with a different number of participants, and that can 
easily manage membrane compartments. In particular, we can encode the 
MBD version of Brane calculi with the link-calculus, and we can also offer 
similar encodings for further versions of the calculus, based on different sets 
of primitives. Our flat calculus puts its set of low-level primitives at our 
disposal to simply implement the various more abstract primitives of Brane. 
We have an operational correspondence result that shows the correctness 
of the proposed encoding. In [8], Cardelli shortly addresses the problem 
of encoding Brane Calculi in BioAmbients [25], by observing some critical 
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aspects such as atomicity of actions and the difficulty in managing the 
interactions through the ambient structure. The encodings presented here 
make us confident that such limitations do not apply to a flat and multiparty 
calculus like the link-calculus. 

Furthermore, we can directly model other biological interactions, like 
those involved in the Receptor-mediated Endocytosis, without passing from 
another process calculus. The simplicity and expressivity of our encoding of 
this mechanism encourages us to think that link-calculus can be adopted 
to provide straight models of biological processes. In future work, we would 
like to enrich the link-calculus in order to take into account also the shape 
and geometry of the biological elements, that influence their interactions. 

It would be interesting to further investigate and test the encoding 
capability of link-calculus, by considering other compartment-based calculi, 
such as BioAmbients [25], the calculus for proteins and cells [21], introduced 
as an extension of «-calculus [15], Bio-PEPA [14], the Calculus of Looping 
Sequences [2], Beta Binders [23], the Shape calculus [3]; and Membrane 
Systems [10, 11]. 
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